Understanding surface-adsorption of proteins: the Vroman effect 
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It is now well accepted that cellular responses to materials in a biological medium reflect greatly 
the adsorbed biomolecular layer, rather than the material itself. Here, we study by molecular 
dynamic simulations the competitive protein adsorption on a surface (Vroman-like effect), i.e. the 
non-monotonic behavior of the amount of protein adsorbed on a surface in contact with plasma as 
a function of contact time and plasma concentration. We show how the effect can be understood, 
controlled and inverted. 

PACS numbers: 87.10.Tf, 81.16.Fg, 87.85. J- 



When nanoparticles are in contact with blood plasma, 
or other biological fluids, biomolecules rapidly coat the 
bare surface in a relatively selective manner [l| . It is in- 
creasingly accepted that the early biological responses 
to nanoparticles will be determined by the adsorbed 
biomolecules rather than the pristine surface alone 0, Q . 
Because of their size 0, 0] nanoparticles are trafficked by 
active transport processes throughout the organism, us- 
ing the information from the protein sequences associated 
with the surface of nanoparticles. Unlike the situation 
for flat macroscopic surfaces say of medical implants, for 
nanoparticles the protein environment changes in differ- 
ent compartments of cells and organs, as the nanoparticle 
traffics. This has lent urgency to the modern interest in 
understanding the phenomenon at a more fundamental 
level [H. Still, we can learn a lot from an understanding 
of the process for flat surfaces [1] . Studying the adsorp- 
tion of Fibrinogen on a surface in contact with blood 
plasma, Vroman found that the surface concentration of 
Fibrinogen shows a maximum at an intermediate contact 
time, indicating that Fibrinogen is replaced with time by 
one or more families of different proteins Q. The phe- 
nomenon is not specific to Fibrinogen, but is a general 
effect for many other proteins 0- The plasma proteins 
compete for the occupation of the surface, resulting in a 
sequential competitive adsorption, known as the Vroman 
effect. 

The effects depends on numerous factors such as the 
plasma dilution, the temperature, and the specific surface 
chemistry . In highly concentrated plasma, the sequen- 
tial adsorption takes place in seconds, but it takes sev- 
eral minutes when the plasma concentrations has phys- 
iological values 0. The more hydrophobic the surface, 
the stronger is the protein adsorption, eventually induc- 
ing irreversible adsorption Here we study the ef- 
fect by Molecular Dynamics (MD) of a model protein 
solution in contact with glass. We consider the three 
most abundant proteins in human blood: Albumin (Alb), 
Immunoglobulin-7 (IgG) and Fibrinogen (Fib) , for which 
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FIG. 1. Schematic representations of different proteins ad- 
sorbed on the surface, (a) Alb is a globular protein charac- 
terized by the radius Ra- In all the panels the continuous 
line represents the surface profile. IgG or Fib are represented 
as ellipsoids that interact with each other along the long axis 
Rs,i (b), or along the short axis Ri (c). In the last case they 
are partially deformed (denatured) by the surface adsorption, 
represented here by a partial overlap of the ellipsoids with the 
surface, (d) At a random moment, the adsorbed ellipsoids can 
have different arrangements. 



competitive adsorption has been observed . We show 
that the sequence of adsorption can be explained in terms 
of the relative bulk concentrations, difusivities and sur- 
face affinities of the proteins and that by thermal or 
chemical energy-depletion is possible to control and in- 
vert the effect. 

Alb is a globular protein, with an almost spherical 
shape. We model Alb- Alb interaction as 
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where Ra is the radius of Alb and r the protein-protein 
distance (Fig. [TJi) . Attraction among protein is not in- 
cluded at this level of description, as it is small compared 
to protein-surface interaction Despite this rough 

zero-order approximation, our results support a posteri- 
ori this assumption within the approximations of our ap- 
proach. IgG and Fib in their folded conformation have 
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non-spherical shapes. In particular, the IgG structure 
resembles a 7 and the Fib resembles an elongated ellip- 
soid. They both can be approximated by ellipsoids with 
two principal axes of rotation along which they interact 
with each other (Fig.[T}D-d). This is encoded through the 
protein-protein potential, within the same protein family, 
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where the index i = I,F stands for IgG (/) and Fib (F), 
with a hard-core radius Ri and a soft-core radius Rs,i 
13|. Interaction between pairs of proteins of different 
families are given by Eq.([2]) with parameters Ri and Rs.i 
equal to the averages of the corresponding parameters for 
each family, and with Ra — Rs.A for Alb. 
The protein-surface interaction is given by 
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where is the attractive energy between the surface and 
a protein of the family i, and = i?i/2^/^, with i = 
A, /, F, is the maximum approach distance between each 
protein and the surface. 

For each family of proteins i, we set the soft-core ra- 
dius Rs,i = Rh the hydrodynamic radius, determined 
experimentally from the diffusion coefficient D through 
the Einstein- Stokes equation D = g^^^, where 77 is the 
viscosity of the medium, under the assumption that the 
proteins can be approximated by a sphere. The hard- 
core radiuses Ri are set by imposing for each protein 
that the experimental surface concentration corresponds 
to the close packing configuration [llj. These conditions 
give Ra = 3.55 nm, Ri = 4.9 nm, Rp = 6.58 nm 
Rs,i = 5.51 nm and Rs.f = 11 um. Protein masses 
Ma = 67 KDa, Mi ^ 150 KDa, Mp ^ 340 KDa, neces- 
sary to determine the time scales, arc known from exper- 
iments [loj . Protein-surface attraction energy can be 
calculated from the adsorption rate constants [ll| . These 
rates are proportional to the probability for a protein i 
to attach to the nearby surface 
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However, the in physical units are not known a priori. 
Hence, we consider the relative probabilities for differ- 
ent proteins ^ oc exp ^^^^^ , from which is possible to 
determine the values of the different energies as 
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adopting ea for Alb as the energy units. We set €a by 
comparing our simulations results with experiments at 
ambient temperature, and get e/ = 2.79 ca and ef = 



6.08 EA- We express all the results in terms of the Alb 
units: M* = jj^ for the mass, R* = for lengths 

and e* = for the energies, t* = ( ^^^^^A ' for the 

time, T* = for the temperature, and AE/eA for 

the energy changes of the buffer. In experiments AE is 
controlled by adding sodium azide, or other depletion- 
energy chemical agents, to the protein solution [12|. In 
the following we drop the * for sake of clarity. 

We perform MD simulations at constant T, constant 
volume V and constant number of proteins Ni, in a par- 
allelepiped with two square faces and four rectangular 
faces. A square face is occupied by the attractive sur- 
face, the other by a wall interacting with the proteins 
through the repulsive part of the ^24, 12 potential. We 
apply periodic boundary conditions (pbc) along the four 
rectangular faces. The volume concentrations of proteins 
is taken to match the average concentrations of the hu- 
man plasma, with ca ~ 4.25 g/dl, c/ = 1.25 g/dl and 
cp = 0.325 g/dl, at Xp ~ 100% plasma concentration 
in blood. When a protein is adsorbed on (released by) 
the surface, we keep its bulk concentrations constant by 
inserting (deleting) a protein of the same family in a 
randomly-chosen empty (occupied) space of the box. 

Experiments are usually carried out for highly diluted 
plasma, at concentration as small as Xp = 0.1%, to slow 
down the adsorption rate to minutes or hours, allowing 
precise measurements. However, such low rates would 
decrease the statistics of our MD simulations. We, there- 
fore, perform our simulations in conditions that are closer 
to those of practical interest, with Xp as high as 100%, 
50% and 25%, by considering different sizes of the sim- 
ulation box while keeping constant the initial number of 
proteins, their relative proportions, and the size of the 
adsorption surface. 

For each Xp we average the results over fourteen in- 
dependent runs, starting from independent initial config- 
urations that have been equilibrated by applying pbc in 
any direction. We find that protein surface concentra- 
tions are non-monotonic in time (Fig. [J). For any 
considered Xp, Alb is the first protein that reaches the 
surface, inducing an increase of C^. When the second 
fastest and second most affine protein, IgG, diffuses to 
the surface, it displaces Alb, leading to a decrease of Ca 
and an increase of Cf . Finally Fib, which is the slowest 
and most affine protein to the surface, takes over decreas- 
ing cf and increasing Cp. Each Cf saturates toward 
an equilibrium value at long times, while the total sur- 
face concentration of proteins is saturated at early times. 
This behavior qualitatively reproduces the Vroman ef- 
fect, apart from the behavior of Fib that here is mono- 
tonic, while in experiments has a maximum due to the 
competitive adsorption with heavier and more surface- 
affine plasma proteins not included in our model [l^ . 
The only effect of reducing Xp is a slowing down in the 
dynamics of the process, as observed in experiments 15 1. 
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FIG. 2. Simulations at T = 300 K and (a) Xp = 100% and (b) 
Xp = 25% show that surface concentration of Alb (O), 
Cf of IgG (□) and Cp of Fib (A) are non-monotonic with 
time, while their sum is (V). At Xp — 50% (not shown) we 
find the same behavior with time-scales intermediate between 
those in (a) and (b). Bulk concentrations are as indicated in 
the text. Errors are smaller than symbol sizes. 



The model allows us to understand the sequence of 
surface occupation as a consequence of the competition 
between the smaller, but less affine, proteins with the 
more affine, but bigger, proteins. For example, we test 
that by increasing the Alb affinity, or artificially setting 
to the same value all the diffusion constants, the effect 
disappears. Therefore, affinity and hydrodynamic radius 
are the relevant protein parameters for the effect. 

Next, we study how energy depiction of the protein so- 
lution affects the sequence of adsorption. Here, for sake of 
simplicity, we decrease T, reducing the kinetic energy of 
the solution, but neglecting possible effects of the protein 
stability. We find (i) that, although the surface affinity 
of Fib is stronger than that for IgG, the latter becomes 
the dominant protein adsorbed on the surface for long 
time scales; (ii) that, by changing Xp, the time scale 
of the process becomes longer, but the inversion of the 
protein concentration is always present (Fig. [3]). Hence, 
the energy depletion leads to an inversion of the Vroman 
effect. 

By comparing the results at different energies, fesT, 
and same Xp (Fig. [2E]) , we observe only a week energy- 
dependence of the times at which each Cf reaches its 
maximum. Hence, the time-scales of the process are 
mainly controlled by the total plasma concentration Xp. 

Once we have understood that the protein layer cov- 
ering the surface is controlled by the energy depletion of 
the system, it is interesting to ask if a sudden change 
of external conditions could induce a different composi- 
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FIG. 3. Surface concentrations Cf as function of time for 
T = 120 K at (a) Xp = 100% and (b) Xp = 25%. At 
long time, Cf > Cp, with an inversion with respect to the 
standard conditions in Fig. [2] where C| > Cf . We find the 
same qualitative behavior at Xp = 50%, not shown. Errors 
and symbols are as in Fig. (2] 
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FIG. 4. The surface concentrations Cf , as a function of time 
for Xp = 100%, is drastically affected when the system un- 
dergoes a sudden change from en energy-depleted condition 
to a normal condition. The vertical dashed line marks the 
time to of the change. We control the energy of the solution 
by changing the external parameter T from T = 120 K to 
T = 300 K. Errors and symbols are as in Fig. [2] 



tion of this layer, determining different biomimetics sur- 
face properties. This situation could occur, for example, 
when a medical device is manipulated in a bioenviron- 
ment whose composition is externally controlled during a 
surgery [isj . In particular, we study the case in which the 
system is first equilibrated under energy-depleted condi- 
tions and subsequently undergoes a sudden change that 
reestablishes the normal conditions (Fig. |4]). 

At short times the energy-depleted system evolves until 
the equilibrium concentrations are reached. Under these 
conditions, as discussed (Fig. [3]), the dominant protein 
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is IgG instead of Fib. At time to we switch to normal 
conditions, forcing the system out of equihbrium. As a 
consequence, the system re-enters a transitory situation 
in which the concentrations Cf evolve until they reach 
their new equilibrium values at long times. In the spe- 
cific case considered here, we observe a fast change in the 
surface concentrations, with Cp of Fib overcoming Cf 
of IgG, being the first, under normal conditions, more 
stable on the surface than the second. The final equi- 
librium concentrations are reached at large times. We 
observe also a sudden change in of Alb, between the 
two equilibrium concentrations characteristics of the two 
values of the external parameters T. However, always 
equilibrates to a value that is smaller then Cf and Cp, 
consistent with its long-time values in Fig. [JISl By de- 
creasing Xp, we find the same qualitative behavior for a 
sudden energy- change, but with the transient regime ex- 
tending to longer times, consistent with Fig. [3] Hence, at 
experimental values of Xp the switching behavior would 
occur on time scales that are comparable to those char- 
acteristic of the Vroman effect. 

We remark that our predictions about inverting the 
Vroman effect by changing the experimental control pa- 
rameters should hold only if the protein adsorption on 
the surface is reversible. If the adsorption is, instead, ir- 
reversible the change of external parameters should not 
lead to a new composition of the protein layer. Indeed, 
under many practical conditions of interest for blood 
plasma, it would appear that the binding is indeed mostly 
irreversible 0, [3] ■ Hence, the switching protocol pro- 
posed above represents a possible experimental way to 
evaluate how strongly irreversible is the adsorption pro- 
cess on a specific surface. 

For an irreversible adsorption process, our findings pre- 
dict that by appropriately controlling the parameters of 
the protein solution, such as the amount of depleted en- 
ergy, it is possible to engineering a specific biomimetic 
covering of a surface. Due to the irreversibility, the pro- 
teins, once adsorbed, cannot desorb from the surface, 
even if the external conditions are modified. Therefore, 
it is feasible to cover a device surface with any desired 
protein composition, targeted to a specific biomimetic 
property, by selecting an appropriate initial condition. 
Subsequently, the device could be used under physiolog- 
ical conditions with no further changes of the protein 
cover and its biomimetic properties. 

In conclusion, we study, by MD simulations of a coarse- 
grained model, the Vroman effect for three kinds of pro- 
teins: Alb, IgG and Fib. We show that the effect is the 
consequence of the differences among the proteins proper- 
ties (mass, size and affinity) which regulate their diffusive 
behavior and their interaction with the surface. These 
differences lead to a process of competitive adsorption 



on a surface, in which the different families of proteins 
occupy sequentially the surface, replacing each other, un- 
til an equilibrium situation is reached. By decreasing the 
total concentration of protein in the solution, keeping the 
relative concentrations fixed, the time scales of the pro- 
cess increase and the maxima of surface concentration for 
each family of proteins occur at longer times. 

We find that the protein surface concentrations at equi- 
librium depend on external control parameters. In par- 
ticular, we find that energy depletion induces a drastic 
change in the composition of the covering protein-layer, 
leading to an inversion of the Vroman effect. Our re- 
sults show that the inversion can be used to quantify 
how strongly irreversible is the process of surface adsorp- 
tion of the proteins, an information useful in studies of 
thromboembolic events 171. Furthermore, these results 



suggest the possibility of engineering the composition of 
the protein layer covering a surface in a controlled way, 
a feature particular relevant in biomimetic applications. 
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